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Electron spin relaxation in bulk III-V semiconductors is investigated from a fully microscopic ki- 
netic spin Bloch equation approach where all relevant scatterings, such as, the electron-nonmagnetic- 
impurity, electron-phonon, electron-electron, electron-hole, and electron-hole exchange (the Bir- 
Aronov-Pikus mechanism) scatterings are explicitly included. The Elliot- Yafet mechanism is also 
fully incorporated. This approach offers a way toward thorough understanding of electron spin re- 
laxation both near and far away from the equilibrium in the metallic regime. The dependences of 
the spin relaxation time on electron density, temperature, initial spin polarization, photo-excitation 
density, and hole density are studied thoroughly with the underlying physics analyzed. We find that 
these dependences are usually qualitatively different in the non-degenerate and degenerate regimes. 
In contrast to the previous investigations in the literature, we find that: (i) In n-type materials, the 
Elliot- Yafet mechanism is less important than the D'yakonov-Perel' mechanism, even for the narrow 
band-gap semiconductors such as InSb and InAs. (ii) The density dependence of the spin relaxation 
time is nonmonotonic and we predict a peak in the metallic regime in both n-type and intrinsic 
materials, (iii) In intrinsic materials, the Bir-Aronov-Pikus mechanism is found to be negligible 
compared with the D'yakonov-Perel' mechanism. We also predict a peak in the temperature de- 
pendence of spin relaxation time which is due to the nonmonotonic temperature dependence of the 
electron-electron Coulomb scattering in intrinsic materials with small initial spin polarization, (iv) 
In p-type III-V semiconductors, the Bir-Aronov-Pikus mechanism dominates spin relaxation in the 
low temperature regime only when the photo-excitation density is low. When the photo-excitation 
density is high, the Bir-Aronov-Pikus mechanism can be comparable with the D'yakonov-Perel' 
mechanism only in the moderate temperature regime roughly around the Fermi temperature of elec- 
trons, whereas for higher or lower temperature it is unimportant. The relative importance of the 
Bir-Aronov-Pikus mechanism decreases with the photo-excitation density and eventually becomes 
negligible at sufficiently high photo-excitation density. The effect of electric field on spin relaxation 
in n-type III-V semiconductors is also studied with behaviors very different from those in the two- 
dimensional case reported. Finally, we find good agreement of our calculation with the experimental 
results. 



PACS numbers: 72.25.Rb, 71.70.Ej, 71.10.-w, 72.20.Ht 
I. INTRODUCTION 



Semiconductor spintronics, which aims at utilizing or 
incorporating the spin degree of freedom in electron- 
ics, has attracted much interest jii^i^ During the last 
decade, the fast developments of techniques of coher- 
ent manipulation of electron spins via optical or elec- 
trical methodsii^i^iii^i^iiaiiiii^ii^ have intrigued a lot of 
studies>2. Many of the important findings are developed 
in bulk GaAs or GaAs epilayers, where the spin relax- 
ation time (SRT) was found to be as long as 130 ns^i^ 
and the spin diffusion length was reported as large as 
100 /im.~ Remarkably, it has been found that the SRT 
can vary by more than three orders of magnitude with 
temperature or electron density^. The relevant spin re- 
laxation mechanisms for electron system in the metal- 
lic regime have been recognized for a long time as: (i) 
the Elliot- Yafet (EY) mechanism in which electron spins 
have a small chance to flip during each scattering due to 



spin mixing in the conduction band;^** (ii) the D'yakonov- 
Perel' (DP) mechanism in which the electron spins de- 
cay due to their precession around the k-dependent spin- 
orbit fields (inhomogeneous broadeningi^) during the free 
fiight between adjacent scattering eventsji^ (iii) the Bir- 
Aronov-Pikus (BAP) mechanism in which electrons ex- 
change their spins with holesiii The hyperfine interac- 
tion is another mechanism which is usually important 
for spin relaxation of localized electrons and ineffective 
in metallic regime where most of the carriers are in ex- 
tended statesi^ii^^^^ 

Despite decades of study, a detailed theoretical in- 
vestigation from a fully microscopic approach in bulk 
system has not yet been performed. Although Song 
and Kim have investigated electron spin relaxation due 
to all the three relevant mechanisms for various con- 
ditions in both n- and p-type semiconductors,^^ in 
their work they use the analytical expressions based on 
single-particle approach which are only applicable for 
the non-degenerate electron system, and hence make 
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the discussion in low temperature and/or high density 
regime questionable. More importantly, the carrier- 
carrier Coulomb scattering, which has been shown to be 
very important for spin relaxation in two-dimensional 
systems,il^^i25^^i28^-30,3i,32,33,34.35 ^as not yet 

been well studied in bulk system. Also, at finite spin 
polarization, the Coulomb Hartree-Fock (HF) term acts 
as an effective longitudinal magnetic field, which has 
been demonstrated to be able to increase the SRT by 
more than one order of magnitude when the initial spin 
polarization is high in two-dimensional electron system 
(2DES),^i26i37 However, the effect of the Coulomb HF 
term in bulk system has not been investigated. An- 
other issue is that the commonly used analytical for- 
mula for spin relaxation due to the BAP mechanism is 
based on the elastic scattering approximation, which has 
been proved to be invalid for low temperature due to 
pretermission of the Pauli blocking of electrons very re- 
cently by Zhou and Wu.^^ Consequently, the BAP mech- 
anism has been demonstrated to be unimportant in two- 
dimensional system at moderate and high excitation den- 
sity first theoretically^^ and then experimentally,-* which 
is in stark contrast with the common belief in the liter- 
ature. Whether it is still true in bulk system remains 
unchecked. Furthermore, in most works only the short- 
range electron-hole exchange interaction was considered, 
while the long-range part was ignored J^'^^i^'^ All these 
questions suggest that a detailed fully microscopic inves- 
tigation is needed. In this work, we perform such a study 
from the fully microscopic kinetic spin Bloch equation 
(KSBE) approach^^i25^i27^^^^di We focus on 
the metallic regime where most of the carriers are in ex- 
tended states. We restrict ourselves to the zero magnetic 
field case 4^ 

Previously, the KSBE approach has been ap- 
plied extensively to study spin dynamics in semi- 
conductor nanostructures in both Markovian and 
non-Markovian limits and in systems both near 
and far away from the equilibrium (under strong 

static or terahertz electric field or with high spin 
polarization) ^^i2i^>26^>28|29^^^^did3d4 ^he 

KSBE approach has been demonstrated to be success- 
ful in the study of spin relaxation in semiconductor 
quantum wells where good agreements with experiments 
have been achieved and many predictions have been 
confirmed by experiments^'^'^i^'^iiSS^ds^di This 
approach has also been applied in the investigation of 
spin relaxation/dephasing in bulk GaAs many years ago 
by Wu and Ning, where only the electron-nonmagnetic- 
impurity and electron-phonon scatterings are included4^ 
In this work, we include all the scatterings, especially 
the electron-electron Coulomb, electron-hole Coulomb, 
and electron-hole exchange scatterings together with the 
EY mechanism which were not considered in Ref. IZsl . 

An important goal of this work is to find the domi- 
nant mechanism in different parameter regimes and for 
different materials from the fully microscopic KSBE ap- 
proach. Previous investigations in the literature indicate 



that for n-type III-V semiconductors, the spin relaxation 
is mostly dominated by the DP mechanism, except at low 
temperature where the EY mechanism is most important; 
for intrinsic and p-type III-V semiconductors, the BAP 
mechanism dominates at low temperature when the hole 
density is high, whereas the DP mechanism dominates in 
other regimes. ^i^i^i^ In the present work, from the fully 
microscopic KSBE approach, however, we find that the 
EY mechanism is less important than the DP mechanism 
in rt-type III-V semiconductors even for narrow band-gap 
semiconductors, such as InAs and InSb. For p-type III-V 
semiconductors, we find that the BAP mechanism domi- 
nates spin relaxation in the low temperature regime only 
when the photo-excitation density is low enough. How- 
ever, when the photo-excitation density is high, the BAP 
mechanism can be comparable with the DP mechanism 
only in the moderate temperature regime roughly around 
the Fermi temperature of electrons, and for higher or 
lower temperature it is unimportant. The relative impor- 
tance of the BAP mechanism decreases with the photo- 
excitation density and eventually becomes negligible for 
sufhciently high photo-excitation density. For intrinsic 
III-V semiconductors, the BAP mechanism is negligible. 

An important method of spin injection is the hot- 
electron spin injection where high electric field is 
applied?^ Moreover, in 2DES, the spin relaxation can 
be effectively manipulated by the high in-plane electric 
fieldi^i^^iS^ In this work, we also study spin relaxation in 
n-type bulk semiconductors under high electric field. We 
show that there is some essential difference of the electric 
field effect on spin dynamics between 2DES and the bulk 
system. Using GaAs as an example, we demonstrate that 
the electric field dependence of spin lifetime can be non- 
monotonic or monotonic depending on the lattice tem- 
perature and the densities of impurities and electrons. 
The underlying physics is analyzed. The results indi- 
cate that the spin lifetime can be effectively controlled 
by electric field. 

This paper is organized as follows: In Sec. II, we intro- 
duce the KSBEs. In Sec. Ill, we study the spin relaxation 
in n-type III-V semiconductors. In Sees. IV and V, we 
investigate the spin relaxations in intrinsic and p-type 
III-V semiconductors, respectively. We study the effects 
of electric field on spin relaxation in n-type III-V semi- 
conductors in Sec. VI. Finally, we conclude in Sec. VII. 



II. KSBES 

The spin dynamics is studied by solving the micro- 
scopic KSBEs derived via the nonequilibrium Green func- 
tion methodr^iS^i^i^ 

dtPk = 9t/5k|coh + 9t/5k|drift + 9t/5k|scat- (1) 

Here pk is the single particle density matrix with the 
diagonal terms representing the distributions of each spin 
band, and the off-diagonal terms denoting the correlation 
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of the two spin bands. The coherent term is given by 

9tPk|coh = -iMk) ■ I + SHF(k), pk]. (2) 

Here [ , ] is the commutator and r2(k) = r2D(k) + r2s(k) 
with 

nD(k) = 27D (fc, (fc2 _ kl) , ky {kl ~ kl) , fc. {kl - kl)) 

(3) 

due to the Dresselhaus spin-orbit couphng (SOC)^"^ and 
f2s(k) =2/3(fc,,-fc^,0) (4) 

due to the strain induced SOCJ-ii^i^ SHF(k) = 
— X^k' ^-k'Pk' is the Coulomb HF term of the electron- 
electron interaction. Previously, it was found that the 
Coulomb HF term serves as a longitudinal effective mag- 
netic field which increases with the initial spin polariza- 
tion. The effective magnetic field can be as large as 40 T, 
which blocks the inhomogeneous broadening of the k- 
depcndent spin-orbit field and reduces the spin relaxation 
due to the DP mechanism in 2DES.^'^ However, the effect 
of the HF effective magnetic field on spin relaxation in 
bulk system is still unknown, which is one of the goal of 
this paper. The drift term is given by,^^ 

9t/5k|drift = -eE • Vk/5k, (5) 

where e is the electron charge (e < 0) and E is the 
electric field. As the hole spin relaxation is very fast 



(~100 fs) )^4,55 pj^g ^g^T^ assume that the hole system is 
kept in the thermal equilibrium state where the hole 
distribution (/k„J is described by the Fermi distribu- 
tion. The scattering term 9tpk|scat contains the contri- 
butions from the electron-impurity scattering 9t/5k|ei, the 
electron-phonon scattering 9tyOk|ep, the electron-electron 
scattering 9t/5k|cc, the electron-hole Coulomb scatter- 
ing 9tPk|oh, and the electron- hole exchange scattering 

<9tPk|ox, 

9t/5k|scat = 9tPk|oi + 9tPk|op + 9tPk|oo 

+ 9t/5k|eh + 9tPk|ex- (6) 

These terms read 
9t/5k|ei = -7r^n,ZfFk_k''5(£k' -ek)(Ak,k'/5k'Ak',k 

k' 

X - Ak^k'P^, Ak',k/5^ ) + H.c, (7) 

9t/5k|op = X! I^^A,k-k'P<5(±WA,k-k' + Ek' - £k) 

A,±,k' 

X (^tk-k'^k,k'P^,Ak',kP^ - ^Jk-k'^k,k' 

X /5^,Ak',kPk) +H.C., (8) 



9tPk|eo = X! ^k-k'^(£k' -£k + £k" -£k"-k+k') 

k',k" 

- Ak.k'Pk'Ak'.kPkTr ^Ak",k"-k+k'/5k"_k+k'^k"-k+k',k"/5k" 



Ak.k'Pk'Ak'.kPkTr (Ak",k"-k-i-k'Pk"-k-i-k'^k"-k-i-k',k"/5k") 
-H.c, (9) 



k' ,k'' ,m,m' 



^k-k''^(£k' £k + £k"m £k"-k+k'm') ^k.k'Pk'^k'.kPk l^'^+k'm' P/k"-k-|-k 



(l .fu."m) ^k,k'Pk'^k',kPk I'^'^+k'm'P (l .fu."-k+k'm')fk."m + H.C 



(10) 



• E 

k' ,k" ,m,m' ,x— =t 



<^(£k' — £k + £k"m ^ £k"-k+k'm') ^xPk"^-xPk l^k"-k+k'm' l^/k"-k+k'm' (l ./k"m) 



*xPk'*-xPk l'^k"-k+k'm'P (l /k"-k+k'm') /k"m +H.C.. 



(11) 



In these equations, Ak,k' — ^ ~ i^d^ x k') • o- de- 
scribes the spin-mixing due to the conduction-valence 
band mixing which originates from the FY spin re- 
laxation mechanismiiiii Here ~ „ with 

^ ~ A^^+E ■ ^9 ^so are the band-gap and the 



spin-orbit splitting of the valence band, respectivelyji 
is the impurity density. Zi is the charge number of the 
impurity which is taken to he Zi = 1 throughout the 
paperi^ ek = k'^/{2mc) with rric being the conduction 
band effective mass. Vq is the screened Coulomb po- 
tential where the screening is treated within the random 
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phase approximation (RPA)^i^i^ 



§)), 



(0) 



l-Kl"^P(l)(q)^ 



(12) 



where 



E 

k,o- 



/k+q,(T — /] 



k.cr 



£k 



E 1^. 



,771' 



fh 



k.m 



k+q.rn' k/ 



(13) 



Here V^^'' = e^/(eo'«o9^) is the bare Coulomb potential 
with eo and kq representing the vacuum permittivity and 
the static dielectric constant, respectively; /k.o- is the 
electron distribution on cr-spin band; /^m is the hole 
distribution on the hole m-spin band. ejj„j = A:^/(2mJ^) 
stands for the hole energy dispersion with m*^ being 
the hole effective mass. For heavy-hole (to = ±|), 
TO^ — Too/ (71 — 272); while for hght-hole (to = ±|), 
m*^ = Too/ (71 + 272). Here 71 and 72 are the param- 
eters of the Luttinger-Kohn Hamiltonian in the spher- 
ical approximation^^ and toq is the free electron mass. 
Note that in these bands the spins are mixedi^ Con- 
sequently, there are form factors in the electron-hole 
Coulomb scattering: \T^Z-? = where 
|^m(k))'s are the eigen-states of the hole Hamiltonian 
which can be found in Ref. [ssl . M^^q is the matrix el- 
ement of electron-phonon interaction with A being the 
phonon branch index, which is further composed of the 
electron-longitudinal-optical(LO)-phonon and electron- 
acoustic-phonon interactions. The expressions of Af^.q 
can be found in Ref. [isl . ojaq is the phonon energy spec- 

^fq . ^ ■ 

number at lattice temperature. 

For electron-hole exchange scattering, the form factor 

l'^k"-k+k'm' 1^ ^1- d) comes from the electron-hole 
exchange interaction. The Hamiltonian for the electron- 
hole exchange interaction consists of two parts: the 
short-range part fl"sR and the long-range part /Jl^.— 1^ 
The short-range part contributes a term proportional to 
~\. \^d{^\'^ ' ^' where Ai?sR is the exchange splitting of 
the exciton ground state and |03d(O)P — l/(7raQ) with oq 
being the exciton Bohr radius. S and s are the hole and 
electron spin operators respectively. The long-range part 



trum. = Nxq + ^ i 5 with iV^q being the phonon 



gives a term proportional to 



3 AEu 



-(Mil + M,s, 



8 |03d(O)|^' 

^M-S+ + iAf+s_) where Ai?LT is the longitudinal- 
transverse splitting; Mi, Mz, M_, and M+(= m1) are 
operators in hole spin space. s± — Sx i: isy are the 
electron spin ladder operators. The expressions for M's 
can be obtained from the expressions of i?LR in Ref. Isol . 
Specifically, the spin-flip matrix M_ in the spin unmixed 
base can be written as (in the order of 



1 



2Kl 



\KzK_ 



3"- 



-fi^.i^- 



.(14) 



Here K = k' -I- k" is the center-of-mass momentum of 
the interacting electron-hole pair. K± = ± iKy and 
= + Ky. Mil + MzSz corresponds to the spin- 
conserving process which is irrelevant in the study of spin 
relaxation. Summing up the contribution from the short- 
range and long-range parts and keeping only the spin-flip 
terms, one obtains 



(±) k'r 



with = 



(e™(k')|:7(±)|6n'(k'-q)), (15) 



AEr 



16 |</.3D(0)r 



^M+ 



1 A^SR 



4 |03d(O)| 



7 5*4 



^4 



Sx ± iSy are the hole spin ladder operators. It should be 
noted that in GaAs, Ai^LT — 0.08 meV is four times as 
large as AEsr = 0.02 meV,^i^ Thus in GaAs the long- 
range part of electron-hole exchange interaction should 
be more important than the short-range part. Later in 
this paper, we will compare the contributions from the 
short-range and long-range parts to show that the long- 
range part is much more important than the short-range 
part in GaAs. We will also compare the results from the 
fully microscopic KSBE approach with those from the an- 
alytical formula widely used in the literaturefi'^'^'^^iiS 

In summary, all relevant spin relaxation mechanisms 
have been fully incorporated in our KSBE approach. By 
solving the KSBEs [Eq. ([T])], one obtains the temporal 
evolution of spin density matrix /3k. After that, the time 
evolution of macroscopic quantities, such as, the electron 
spin density along the z-axis, Sz(t) = Tr [szPk(t)], are 
obtained. By fitting the decay of Sz with an exponential 
decay, one obtains the SRT. The SRTs under various con- 
ditions are studied with the underlying physics discussed. 
The material parameters used are listed in Table I. The 
parameter for strain-induced SOC (3 is always taken to 
be zero unless otherwise specified. The numerical scheme 
is laid out in Appendix A. The error of our computation 
is less than 5 % according to our test. 

We first compare our calculation with experiments. In 
Appendix B, we compare the SRTs calculated from the 
KSBEs with those from experiments in Refs. 0, [H and 
40. We find good agreement with experimental data in 
almost the whole temperature or photo-excitation den- 
sity range with only one fitting parameter 70-^^ This 
demonstrates that our calculation has achieved quantita- 
tive accuracy in modeling the spin relaxation in metallic 
regime. 

In the rest part of this section, we briefiy comment 
on the merits of the fully microscopic KSBE approach. 
First, for the DP mechanism. Previously, the widely used 
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TABLE I: Material parameters used in the calculation (from 
Ref. [go! unless otherwise specified) 





GaAs 


GaSb 


InAs 


InSb 


E, (eV) 


1.52 


0.8113 


0.414 


0.2355 


Aso (eV) 


0.341 


0.75 


0.38 


0.85 


rric/mo 


0.067 


0.0412 


0.023 


0.0136 


Ko 


12.9 


15.69 


15.15 


16.8 




10.8 


14.44 


12.25 


15.68 


tJLO (meV) 


35.4 


28.95 


27.0 


23.2 


v,i (10^ m/s) 


5.29 


4.01 


4.28 


3.4081 


Vst (10^ m/s) 


2.48 


2.4 


1.83 


2.284 


D (10^ kg/m^) 


5.31 


5.6137 


5.9 


5.7747 


-(eV) 


8.5 


8.3 


5.8 


14.0 


ei4 (V/m) 


1.41x10^ 


9.5x10** 


0.35x10^ 


4.7x10^ 


7D (eV-A^) 


23.9"''' 


168'' 


42.3'' 


389*' 


A£;lt (meV) 


0.08'= 


0.02'* 






ASsR (meV) 


0.02= 


0.024"= 






71 


6.85 


11.8 


19.67 


35.08 




2.5^ 


4.65^ 


8.83-'^ 


16.27^ 


" Ref. 
" Ref.m. 


Ref. 63. " Ref. 61. Ref. 64. 
^ Obtained from Ref. i60j by i(72 + 73) 



analytical formula for the SRT due to the DP mechanism 



Tup 



En 



(16) 



where Q is a numerical constant which varies around 1 
depending on the dominant momentum scattering mech- 
anism, r„ is the momentum relaxation time, and a — 



2jDy^2m^Eg is a dimensionless constant. This formula 
is derived from 



'DP 



((|n(k)|2_f72(k))r,(k)) 



(17) 



[Tp (k) is the momentum relaxation time of the state with 
momentum k] by performing ensemble averaging over the 
Boltzmann distribution.'^ An important fact about this 
formula is that it is derived in the elastic scattering ap- 
proximation, which artificially confines the random-walk- 
like spin precession due to the k-dependent spin-orbit 
field only within the same energy states. However, in 
the genuine case the inelastic electron-phonon scattering 
(especially the electron-LO-phonon scattering) as well 
as the carrier-carrier Coulomb scattering can be more 
important, and the random spin precession (the "in- 
homogeneous broadening" ) should be fully counted for 
the whole k-space, instead only within the same energy 
states. On the other hand, the KSBE approach, which 
solves the spin precession and the momentum scatter- 
ing self-consistently, takes full account of the inhomo- 
geneous broadening and the counter effects of scatter- 
ing on the inhomogeneous broadening. Moreover, the 
electron-electron scattering (although it does not con- 



tribute to the mobility) as well as the electron-hole scat- 
tering, which have been demonstrated to have important 
effects on spin relaxation/dephasing in two-dimensional 
systems^^^^i^e^^S^^^^^^ are always ne- 



glected in previous studies based on Eq. p6)) . On the con- 
trary, the fully microscopic KSBE approach includes both 
the electron-electron and electron-hole Coulomb scatter- 
ings. In this work, we will discuss the effect of the 
electron-electron and electron-hole Coulomb scatterings 
on spin relaxation. 

Second, for the EY mechanism. The formula com- 
monly used in the literature for the SRT due to the EY 
mechanism reads^ 

2 /I /r,\ 2 



1 

TEY 



A 



kBT 

Ea 



1-V/2Y 1 



1-77/3/ Tj 



(18) 



where the numerical factor A is of order 1 depending 
on the dominant scattering mechanism. This formula, 
which is also based on the elastic scattering approxima- 
tion, has similar problems with those of Eq. (jl6p . On the 
contrary, the fully microscopic KSBE approach, which 
incorporates all the EY spin-flip processes in all relevant 
scatterings)^ fully takes into account of the spin relax- 
ation due to the EY mechanism. 

Third, for the BAP mechanism. The commonly used 
formula for the SRT due to the BAP mechanism for non- 
degenerate holes (assuming all holes are free) is^ 



(19) 



1 2 3 (t-k) 
tbap To vb 

where ob is the exciton Bohr radius, 1/tq = 
(37r/64)Ai?gp,/i?B with E^ being the exciton Bohr en- 
ergy, Tih is the hole density, (uk) = {k/rric) is the average 
electron velocity, and wb = I/I'ti/jQ-b) with mji w rric 
being the reduced mass of the electron-hole pair. For 
degenerate holes,.- 

1 3 .(t^w)fcBr 



TBAP 



6 3 Wk^ 

To VB E^ 



with Ep denoting the Fermi energy of holes. All these 
formulae are obtained within the elastic scattering ap- 
proximation. Previously, it has been shown that for two- 
dimensional system the elastic scattering approximation 
is invalid for low temperature regime where the Pauli 
blocking of electrons is important In bulk system, sim- 
ilar problem also exists^^ However, the KSBE approach 
keeps all the electron- hole exchange scattering terms, and 
thus gives correct results. 

Finally, we point out that for excitation far away from 
the equilibrium, such as excitation with high spin po- 
larization, the spin-conserving scatterings are very im- 
portant for redistribution of electrons in each spin band, 
which also affects the spin dynamics. This effect is auto- 
matically kept in our approach, but missing in the ana- 
lytical formulae. 

Although the analytical formulae based on single 
particle approach [Eqs. ([TC)l -(|20 p ] have some disad- 
vantages, for non-degenerate electron system they still 
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give qualitatively good results. However, quantitative 
analysis based on them should be questioned. It has al- 
ready been extensively demonstrated in two-dimensional 
system (both theoreticalljii^^'^s^iHl^^^^^^ds 
and experimentally22^^^i37^did5d6^S^^^) 
that the above single-particle approach, with even the 
recent developments which exert closer examination 
on both the electron distribution and the energy- 
dependent scattering cross sections, are inadequate in 
accounting for the spin relaxationj^^ The same is true 
for bulk systemi^i22ii2^i2ii2^ Nevertheless, when the 
electron-impurity scattering is much stronger than the 
electron-electron (electron-hole) and electron-phonon 
scatterings [i.e., in degenerate electron (hole) regime in 
n-type (p-type) bulk semiconductors with low excita- 
tion density], the above approach for the DP and EY 
spin relaxations and the later developments^ii^Sii^^ are 
applicable, thanks to the fact that the impurity density 
is larger than or equal to the electron (hole) density 
in n-type (p-type) bulk semiconductors (which is quite 
different from the two-dimensional system). On the 
other hand, the fully microscopic KSBE approach takes 
full account of the DP, EY and BAP spin relaxations 
by solving the kinetic equation directly. It is applicable 
in all parameter regimes in the metallic regime. Fur- 
thermore, the fully microscopic KSBE approach can be 
applied to system far away from equilibrium, such as 
under strong electric field where the hot-electron effect is 



above conclusion fails. A revised expression is obtained 
by replacing k^T with the average kinetic energy (ek)n 



prominent 



24.29.30 



or with high spin polarization, and/or 



hot photo-carriers , ^^^se j-j- jg even valid for system out of 
the motional narrowing regime, or in the non-Markovian 
hmit^^ 



III. SPIN RELAXATION IN n-TYPE III-V 
SEMICONDUCTORS 



Comparison of different spin relaxation 
mechanisms 



In rt-type III-V semiconductors, the BAP mechanism is 
ineffective due to the lack of holes. The remaining mech- 
anisms are the DP and EY mechanisms. Previously, it is 
widely accepted in the literature that the EY mechanism 
dominates spin relaxation at low temperature, while the 
DP mechanism is important at high temperaturejii^ii^^ 
Most studies concerned are based on the analytical for- 
mulae [Eqs. p6|) and p8)) ]. which give 



T'EY 

Tup 



(21) 



From the above equation, one arrives at the conclusion 
that the EY mechanism is dominant at low temperature 
and/or when the momentum scattering is strong (such 
as, for heavily doped samples). However, the electron 
system may enter into the degenerate regime at low tem- 
peratures and/or in heavily doped samples, in which the 



tey 2(5 , 2o 
TDP 3A ^ 



(22) 



with 8 = ^ 1-^/2 ) • This is still correct only qual- 
itatively. 

TABLE II; The factor 6 for IILV semiconductors 



GaAs 


GaSb 


InAs 


InSb 


InP 


e (eV) 0.23 


0.12 


3.6x10"'' 


9.2x10"" 


0.27" 



" The SOC parameter 7d is from Ref. i6S . 
other parameters are from Ref. IgoI . 

In this work, we reexamine the problem from the fully 
microscopic KSBE approach. Let us first examine the 
case of InSb. InSb is a narrow band-gap semiconductor 
where the spin relaxation is believed to be dominated by 
the EY mechanism at low temperature previously . '^'^^i^^ 
The factor Q is much smaller than that of GaAs, which 
indicates the importance of the EY mechanism according 
to Eq. ((22)) (A list of the factor 8 for different materials 
is given in Table II). In Fig. [IJa), we plot the ratio of 
the SRT due to the EY mechanism, tey, to that due to 
the DP mechanism, tdp, calculated from the KSBEs, as 
function of temperature for various electron densities. In 
the calculation, Ui = Ue- Remarkably, it is noted that 
the ratio tey/tdp is always larger than 1, and in most 
cases it is even larger than 10, i.e., the spin relaxation 
in n-type InSb is not dominated by the EY mechanism. 
Moreover, the temperature dependence is not monotonic, 
which is different from the intuition given by Eq. (j22p . 
This can be understood as following: The increase of the 
ratio comes from the increase of (sk) which is understood 
easily. The decrease of the ratio comes from the decrease 
of Tp due to the increase of the electron-LO-phonon scat- 
tering with temperature which becomes important for 
T > 80 K in InSb. Therefore, after a crossover regime, 
the ratio eventually decreases with temperature. More- 
over, for large electron density rig > 10^^ cm"'^, the elec- 
tron system is in the degenerate regime and (ek) changes 
slowly with temperature which further facilitates the de- 
crease of the ratio. To elucidate it clearly, we plot the 
SRT due to the EY mechanism and DP mechanism in 
Fig. [TJb) and (c) respectively. It is seen from Fig. [Tfb) 
that the SRT due to the EY mechanism always decreases 
with temperature due to both the increase of (ek) and 
the enhancement of the scattering. From Fig. [T](c), it is 
further noted that the temperature dependence of the 
SRT due to the DP mechanism is, however, different 
for low density and high density regimes: in low density 
regime, the SRT decreases with the temperature, whereas 
in high density regime it increases. The decreases of SRT 
with temperature at low density is consistent with pre- 
vious studies in the literature.— The increases of SRT 
with temperature at high density, however, is because 
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that the electron-LO-phonon scattering increases with 
temperature faster than the inhomogeneous broadening 
^ ((|r2(k)p — ri^(k))) cx (e^) in the degenerate regime. 
Note that the increase of the SRT with temperature at 
such high temperature (T > 100 K) without any mag- 
netic field has neither been observed in experiment nor 
been predicted in theory. Note that the results for high 
density cases are only given at high temperatures due 
to the limitation of our computation power. However, 
according to the above analysis, the ratio in the low tem- 
perature regime should be larger than its value at 300 K 
for high density cases. For low density cases with tem- 
perature lower than our calculation range, both the mo- 
mentum scattering and the inhomogeneous broadening 
change little with temperature as the electron system is 
degenerate. Therefore the ratio tey/^dp changes slowly. 
Consequently, it does not change the conclusion that the 
EY mechanism is less efficient than the DP mechanism 
in n-type InSb.''^ In contrast, based on Eq. ([^ . Song 
and Kim reported that the EY mechanism is more im- 
portant than the DP mechanism for temperature lower 
than 5 

What makes our conclusion different from that in the 
literature is that most of the previous investigations use 
Eqs. (fT6|) and (fT8)) to calculate tdp and tey^'^ However, 
these equations are only applicable for non-degenerate 
electron system. If it is used in the degenerate electron 
system, it exaggerates the relative efficiency of the EY 
mechanism according to Eq. (|2ip . 

We further examine other two III-V semiconductors: 
InAs and GaAs. Figure [5] shows the ratio tey/tdp for 
different electron densities as function of temperature for 
both InAs and GaAs. It is seen that the ratio is al- 
ways larger than 1 for both InAs^^ and GaAs,— which 
is different from previous results in the literature. 
Especially, for GaAs the ratio is larger than 100, which 
indicates that the EY mechanism is irrelevant for rt-type 
GaAs. The temperature dependence of the ratio is simi- 
lar to that in InSb and the underlying physics is also the 
same. 

In summary, we find that the EY mechanism is less 
efficient than the DP mechanism in n-type GaAs, InAs, 
and InSb. According to the data listed in Table II, we 
believe that the same conclusion holds for other n-type 
III-V semiconductors. 
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FIG. 1: (Color online) n-InSb. Ratio of the SRT due to 
the EY mechanism tby to that due to the DP mechanism 
T"DP (a), tey (b) and tdp (c) as function of temperature for 
various electron densities. Tie 
3 X 10 
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^ (curve with •), 
(curve with A), 



cm (curve with □), 10 
3 X 10^^ cm~"^ (curve with V). The corresponding Fermi tem- 
peratures for those densities are Tp = 144, 300, 670, and 
1390 K respectively. Ui = Ue- 



B. DP spin relaxation 

As the BAP mechanism and the EY mechanism are 
unimportant for spin relaxation in n-type III-V semicon- 
ductors, in this subsection we focus on spin relaxation 
due to the DP mechanism. The effect of the electron- 
electron Coulomb scattering on spin relaxation is stud- 
ied. Previously, it was found that the electron-electron 
scattering plays an important role in two-dimensional 
system^^i^i^i^e^iSS^^^ Especially, it becomes the 
dominant scattering mechanism in high mobility sam- 



ples at low temperature, where other scattering mecha- 
nisms are relatively weak.^-i^i^i^ Consequently Tj,(k) 
in Eq. dTT]) should be replaced by T*{\<i) with l/rp*(k) 
including the electron-electron scattering 1/Tp^(k) [later 
we will use the symbols t* and t^'^ to denote the en- 
semble averaged value] J^i^^i^li^ Remarkably, the SRT 
due to the electron-electron scattering has nonmonotonic 
temperature dependence: in the low temperature (de- 
generate) regime the SRT increases with temperature as 
the electron-electron scattering does; in the high tem- 
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FIG. 2: (Color online) n-InAs (solid curves) and n-GaAs 
(dashed curve). Ratio of the SRT due to the EY mecha- 
nism tey to that due to the DP mechanism tdp as function 
of temperature for various electron densities: Ue — 10^^ cm~^ 
(curve with •), 2 x 10^'' cm"^ (curve with □), 10^** cm~^ 
(curve with A) for InAs, and Ue — 10^* cm""* (curve with V) 
for GaAs (note that the value in the figure has been rescaled 
by a factor of 0.02). The Fermi temperature Tp is 85, 629, 
and 1840 K for InAs as well as 29 K for GaAs respectively. 

Hi = He- 

perature (non-degenerate) regime the electron-electron 
scattering decreases with the temperature, so does the 
SRT Jii^J Thus there is a peak Tc in the temperature de- 
pendence of the SRT which is comparable with the Fermi 
temperaturci ^^i^^ This prediction was confirmed by ex- 
periments very recently^ However, in bulk semiconduc- 
tors, the role of electron-electron scattering in spin relax- 
ation is still unclear. Although Glazov and Ivchenko have 
discussed the problem, they only gave an approximate 
expression of the SRT due to the electron-electron scat- 
tering in the non-degenerate regime, while the relative 
importance of the electron-electron scattering compared 
with other scattering mechanisms is not addressedi^ In 
this subsection, we present a close study on the effect of 
the electron-electron scattering on spin relaxation in bulk 
semiconductors. We use GaAs as an example, where the 
behavior is applicable to all III-V semiconductors given 
that the system is in the motional narrowing regime. 

1. Electron density dependence 

We first discuss the electron density dependence of the 
SRT. It is noted that during the variation of electron 
density, the impurity density also varies as = Ue- This 
is different from the situation in 2DES, where the im- 
purity density can be different from the electron density 
due to the modulation doping. In Fig.[3l[a), we plot the 
SRT as function of electron density for T = 40 K. It is 
noted that, remarkably, the density dependence of the 
SRT is nonmonotonic and there is a peak in the r-rie 
curve. Previously, the nonmonotonic density dependence 
of the SRT was observed in low temperature (T < 5 K) 
measurements, where the localized electrons play an im- 




(cm ) 



FIG. 3: (Color online) n-GaAs at T = 40 K. (a) SRT 
r as function of electron density [ui = rie) from full 
calculation (curve with •), from calculation with only the 
electron-electron scattering (curve with □), with only the 
electron- impurity scattering (curve with A), and with only 
the electron-phonon scattering (curve with V); (b) SRT r as 
function of electron density Ue {rn = rie) for the case with 
strain (curve with •: with both the linear and the cubic SOC; 
curve with A: with only the linear SOC) and the case without 
strain (curve with □). 



portant role and the electron system is in the Mott metal- 
insulator transition area i^'^i^i^^'^^i^^ The localized elec- 
trons have different spin relaxation mechanisms and the 
scatterings of the localized electrons and free electrons 
give rise to the nonmonotonic density dependenceJ^i^ 
It is noted that, up till now, there is no report on the 
nonmonotonic density dependence in the metallic regime 
in 71- type bulk III-V semiconductors. It should be further 
pointed out that the nonmonotonic density dependence 
and the appearance of the peak is a universal behavior in 
the metallic regime in n-type bulk III-V semiconductors. 

We also plot the SRTs calculated with only the 
electron-electron scattering (curve with □), with only the 
electron-impurity scattering (curve with A), and with 
only the electron-phonon scattering (curve with V) in 
the figure to elucidate the role of these scatterings in 
spin dynamics. It is noted that the electron- impurity 
and electron-electron scatterings are the relevant scatter- 
ing mechanisms, while the electron-phonon scattering is 
much weaker than the two as the temperature is low. In- 
terestingly, both the electron-electron scattering and the 
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electron-impurity scattering lead to nonmonotonic den- 
sity dependence of SRT. 

For the SRT due to the electron-electron scattering, 
the nonmonotonic behavior comes from the nonmono- 
tonic density dependence of the electron-electron scat- 
tering. The density and temperature dependences of 
the electron-electron scattering have been investigated in 
spin-unrelated problems. From these works, after some 
approximation, the asymptotic density and temperature 
dependences of the electron-electron scattering time t^'^ 
in the degenerate and non-degenerate regimes are given 

Tp^'^ (X 71^/3 for r<TF, (23) 
T^" (xT^/'^/n^ for r>TF. (24) 

From above equations, one notices that the electron- 
electron scattering in the non-degenerate and degener- 
ate regimes has different density and temperature depen- 
dence. In the non-degenerate (low density) regime the 
electron-electron scattering increases with electron den- 
sity, while the inhomogeneous broadening ~ ((|r2(k)p — 
ri^(k))) changes slowly as the distribution function is 
close to the Boltzmann distribution. The SRT thus in- 
creases with density. In degenerate (high density) regime, 
both Tp'^ and the inhomogeneous broadening increases 
with electron density. Thus, the SRT, which can be es- 
timated as r - l/[(|JI(k)|2 - J72(k))Tp^^]^i2i decreases 
with density. 

For the SRT due to the electron-impurity scattering, 
the scenario is similar: In the non-degenerate regime, 
the distribution function is close to the Boltzmann dis- 
tribution. The inhomogeneous broadening hence changes 
slowly with density. The electron-impurity scattering in- 
creases as (X riiiy^), which increases with the impurity 
density because ly^) changes little with the density when 
the distribution is close to the Boltzmann distribution. 
The SRT hence increases with density. In the degener- 
ate regime, the inhomogeneous broadening increases as 
- ((|f^(k)|2 - 02(k))) (X fc| cx nl. On the other hand, 
the electron-impurity scattering decreases with density 
because it is proportional to ~ niV^^ cx Ue/k^ oc n7'^^^- 
Consequently, the SRT decreases with density. 

Therefore, both the electron-electron and electron- 
impurity scatterings contribute to the nonmonotonic 
density dependence of the SRT. The peak density ric ap- 
pears in the crossover regime, where the corresponding 
Fermi temperature is comparable with the lattice tem- 
perature. A careful calculation gives the peak density 
as ric = 1.4 X 10^^ cm~'^, with the corresponding Fermi 
temperature being 37 K, close to the lattice temperature 
of 40 K. 

We further discuss the electron density dependence of 
the SRT with strain-induced SOC. In Fig.El^b), we plot 
the density dependence of the SRT for T = 40 K un- 
der strain with = rig. We choose (3 = 2.6 meV-Aji^ 
For this value, at low density the SOC is dominated by 
the linear term due to strain. However, in high den- 



sity regime (rig > 2 x 10"'^^ cm""^), the cubic Dressel- 
haus term can surpass the linear term. It is noted that 
the density dependence is also nonmonotonic and ex- 
hibits a peak. The underlying physics is similar: In the 
non-degenerate regime, the SRT increases with the den- 
sity because both electron-electron and electron-impurity 
scattering increase with density and the inhomogeneous 
broadening changes little. In the degenerate regime, the 
inhomogeneous broadening increases as oc /cp oc rie^^, 
whereas both the electron-electron and electron-impurity 
scatterings decrease with the electron density, which thus 
leads to the decrease of the SRT with density. It is fur- 
ther noted that the peak is at ric = 2 x 10^^ cm^^ which 
is larger than the peak density in the case without strain. 
This is because that the inhomogeneous broadening here 
increases as (fc^), which is much slower than that of (fc^) 
in the strain-free case. Nevertheless, the increase of the 
scattering with density remains the same, therefore the 
peak shows up at a larger electron density. This is also 
confirmed in the figure that the SRT with only the linear 
term decreases slower than the one with only the cubic 
Dresselhaus term in the high density regime. However, 
in the low density regime, the SRT in the case with strain 
increases as fast as the one in the strain- free case. This 
is because here the increase of the SRT is due to the 
increase of the scattering, whereas the inhomogeneous 
broadening changes little. 

2. Temperature dependence 

We now study the temperature dependence of the 
SRT. In Fig. [H the SRT as function of temperature is 
plotted for — 10^^ cm^"^. From the figure, it is 
seen that the SRT decreases with temperature mono- 
tonically, which coincides with previous experimental 
results This trend is also the same as that 
in the 2DES with high impurity density.^^ However, for 
high mobility 2DES, which can be achieved by modula- 
tion doping, the temperature dependence of the SRT is 
nonmonotonic and there is a peak Tc around the Fermi 
temperature due to the electron-electron scattering»2I In 
n-type bulk materials, as the impurity density is always 
equal to or larger than the electron density, this peak 
disappears. Nevertheless, similar effect can be obtained 
if one artificially reduces the impurity density. For ex- 
ample when rii — O.Olng, it is seen from Fig. d] that the 
SRT shows nonmonotonic behavior with a peak around 
40 K. 



3. Initial spin polarization dependence: the effect of the 
Coulomb HF term 

We now turn to the dependence of initial spin polar- 
ization on spin relaxation. Previously, it was discovered 
that at finite spin polarization, the Coulomb HF term 
serves as an effective magnetic field along the direction 
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FIG. 4: (Color online) n-GaAs. SRT r as function of T for 
ni — Tie (solid curve with •) and n,; = O.Olrie (dashed curve 
with □). Tie = 10" cm"^. 

of the spin polarization.^"^ The effective magnetic field 
can be as large as 40 T at high spin polarization in 2DES 
which suppresses the DP spin relaxation. This effect 
was first predicted by Weng and W\i^^ and then con- 
firmed by experiments very recently. "^^'^^'^^ However, the 
effect of the Coulomb HF term on spin relaxation in bulk 
system still needs to be evaluated. Here, we present such 
an investigation. 
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FIG. 5; (Color online) n-GaAs. Dependence of SRT r on 
initial spin polarization P for rn — Ue with (without) the 
Coulomb HF term [curve with • (□)] and for rii — O.Olrie 
with (without) the Coulomb HF term [curve with A (V)]. 
ne = 10" cm-=* and T = 20 K. 

The Coulomb HF term for spin polarization along, e.g., 
the z-direction can be written as 

EHF(k) = - ^ Vi_k. (/k'T - /k'i) s.. (25) 

k' 

The corresponding effective magnetic field is along the 
z-axis, 

i?HF(k) = - ^ Vk_k' (/k'T - /k'i) /ffMB. (26) 

k' 

Under this effective magnetic field, the spin preces- 
sion IS blocked and the SRT is elongatedi^i^ Re- 
cently, this effective magnetic field has been probed 



experimentally,- and its effect on spin accumulation 
in 2DES was also discussed theoreticallyi^ 

In Fig. m we plot the SRT as function of the initial 
spin polarization P for rig = 10^^ cm""^ and T = 20 K. 
It is seen that the SRT increases with the initial spin 
polarization. To elucidate the effect of the Coulomb HF 
term, we also plot the results from the calculation with- 
out the Coulomb HF term. The results indicate that the 
increase of the SRT with initial spin polarization is due 
to the Coulomb HF term. However, the increment is less 
than 50 %. Previously, it was shown in 2DES that the 
SRT can increase over 30 times for low impurity case at 
120 K, while less than 3 times for high impurity or high 
temperature case where the scattering is strong.— The 
results can be understood as follows: The SRT under the 
HF effective magnetic field can be estimated as^ 

Ts(P) = Ts{P - 0)[1 + (gMBBHFT;)^] (27) 

where _Bhf is the averaged effective magnetic field. Thus 
the effect of the HF effective magnetic field increases with 
Tp, i.e., the effect is more pronounced for weak scattering 
case, such as the low impurity density case. However, in 
bulk system the impurity density is always equal to or 
larger than the electron density rij > rie, therefore the 
effect of the Coulomb HF term is suppressed. This can 
be seen from the calculation with rii — O.Olng (the artifi- 
cial case). The results are also plotted in the figure. One 
finds that the Coulomb HF term effectively enhances the 
SRT at low impurity density by 40 times. Therefore, due 
to the large impurity density (jii > rie), the SRT is in- 
sensitive to the initial spin polarization in n-type bulk 
III-V semiconductors. We have checked that the con- 
clusion holds for other cases with different temperatures, 
electron densities and materials. 



IV. SPIN RELAXATION IN INTRINSIC III-V 
SEMICONDUCTORS 

In this section we study spin relaxation in intrinsic 
III-V semiconductors. For intrinsic semiconductors, the 
electrons and holes are created by optical excitation, and 
their numbers are equal. Compared to n-type semicon- 
ductors, there are two new scattering mechanisms: the 
electron-hole Coulomb and electron-hole exchange scat- 
terings, where the latter corresponds to the BAP mech- 
anism. Another important property of intrinsic semicon- 
ductors is that the impurity density is very low (we take 
rii = 0), which offers a good platform for demonstrating 
the effect of the many-body carrier-carrier scattering on 
spin relaxation. Moreover, the effect of the Coulomb HF 
term would be enhanced as the electron-impurity scat- 
tering can be eliminated. Our first goal is to compare 
the relative efficiency of the DP and BAP mechanisms. 
After that we also study the temperature and photo- 
excitation density N^x dependences of the SRT. The role 
of electron-hole Coulomb scattering as well as the effect 
of the Coulomb HF term are also addressed. We further 
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FIG. 6: (Color online) Intrinsic GaAs with iVex = 10^^ cm"^. 
(a) SRT T due to the BAP and DP mechanisms as function of 
temperature, (b) SRT due to the BAP mechanism calculated 
from Eq. ([TgJ (dotted curve with A), from the KSBE ap- 
proach with both long-range and short-range exchange scat- 
terings (solid curve with •) as well as from the KSBE ap- 
proach with only the short-range exchange scattering (solid 
curve with □). (c) SRT due to the DP mechanism from full 
calculation (solid curve with •), from the calculation with- 
out the Coulomb HF term (dashed curve with □), and from 
the calculation without the electron-hole Coulomb scattering 
(dotted curve with A). For electrons Tp — 136 K, and for 
holes Tp = 16 K. 



compare the results from the KSBEs with those from the 
widely used analytical formulae [Eqs. ([TO]) and (^0]) ]. The 
initial spin polarization is chosen to be 50 % which cor- 
responds to circularly polarized optical excitation. We 
focus on GaAs, while the situation is similar for other 
III-V semiconductors!^ 



A. Temperature dependence 

In Fig. [6ta), we plot the SRTs due to the DP and 
BAP mechanisms as function of temperature for TVox = 
10^'' cm~^. The SRT due to the BAP mechanism alone 
is calculated by removing the spin precession due to the 
SOC, but keeping all the scattering terms. It is noted 
that the SRT due to the BAP mechanism is larger than 
that due to the DP mechanism by more than one order of 
magnitude, which indicates that the BAP mechanism is 
negligible for intrinsic GaAsi^ Moreover, the spin relax- 
ation due to the DP mechanism increases with tempera- 
ture more rapidly than that due to the BAP mechanism 
at high temperature. This is because the increase of spin 
relaxation due to the DP mechanism mainly comes from 
the increase of the inhomogeneous broadening which is 
proportional to in high temperature (non-degenerate) 
regime. Meanwhile, according to Eq. p9)) . the increase of 
spin relaxation due to the BAP mechanism mainly comes 
from the increase of (vk) which is proportional to T*''^ in 
that regime. 

For a close examination of the BAP mechanism, we 
also plot the SRT limited by the the short-range electron- 
hole exchange scattering calculated from the KSBEs in 
Fig. [DJb) . It is seen that the SRT limited by the short- 
range electron-hole exchange scattering is much larger 
than that limited by both long- and short-range exchange 
scattering. This confirms that the long-range scattering 
is more important than the short-range one in GaAs as 
AEuT is four times larger than Ai?sR. Therefore, pre- 
vious investigation a^'^^1'^^ with only the short-range ex- 
change scattering included are questionable. Moreover, 
to check the validity of the widely used elastic scatter- 
ing approximation, we also compare the results from the 
KSBEs with those from the elastic scattering approxima- 
tion. Under the elastic scattering approximation,^ 
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The SRT is then obtained by averaging over the elec- 
tron distribution. Eqs. p9)) and (pO)) are derived from 
the above equation under some approximations. For ex- 
ample, by including only the short-range exchange scat- 
tering and ignoring the light-hole contribution as well as 
the term of (1 — /k'+qm)' is obtained. To show 

that the elastic scattering approximation fails in the de- 
generate regime, we compare our results with the results 
from Eq. (|28p . For simplicity, we include only the short- 
range exchange scattering. In Fig. [6l[b), we plot the SRT 
obtained from Eq. ([^5]) as the dotted curve. It is seen 
that the result from Eq. agrees well with our result 
from the KSBEs at high temperature, but deviates at low 
temperature. The deviation is due to the Pauli blocking 
of electrons in the degenerate regime, which is neglected 
in the elastic scattering approximation (Tp = 136 K).^" 
We also discuss the effect of the electron-hole Coulomb 
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FIG. 7: (Color online) Intrinsic GaAs with TVcx — 2 x 
10^^ cm~^. SRT r as function of temperature for P — 2 % 
(curve with □) and the SRT as function of initial spin polar- 
ization P for r = 20 K (curve with •) (note that the scale of 
P is on the top of the frame). 



semiconductors, as the impurity density is very low, the 
peak may appear. Indeed, we find that the SRT has a 
peak at Tc ~ 100 K. The peak temperature is compa- 
rable with the Fermi temperature (Tp — 216 K) [Actu- 
ally, our calculation indicates that the peak temperature 
Tc is around Tp/S and lies in the range of (Tf/4,Tf/2) 
depending on the carrier density.] Nevertheless, at high 
spin polarization, such as P = 50 %, the peak disappears 
as indicated in Fig. [SJc). This peak can be observed 
within current technology of optical orientation. How- 
ever, up till now, no such experimental investigation has 
been performed. In Fig. [71 we also plot the SRT as func- 
tion of initial spin polarization at T = 20 K. It is seen 
that the SRT is elongated by 9 times when P is tuned 
from 2 % to 50 %. Therefore, the effect of the Coulomb 
HF term can also be observed in intrinsic materials and 
is more pronounced compared to the n-type case. 



scattering and the Coulomb HF term on spin relaxation 
due to the DP mechanism. In Fig.lHI^c), we plot the SRTs 
due to the DP mechanism obtained from the full cal- 
culation, from the calculation without the electron-hole 
Coulomb scattering, and from the calculation without 
the Coulomb HF term. Let us first examine the effect 
of the Coulomb HF term on spin relaxation. It is seen 
that the Coulomb HF term has important effect on spin 
relaxation only for low temperature case^i (T < 60 K) 
which is consistent with the results in 2DESi^ We then 
turn to the effect of the electron-hole Coulomb scatter- 
ing. It is seen that without the electron-hole Coulomb 
scattering the SRT is larger for T < 60 K but smaller 
for T > 60 K compared with that from the full calcu- 
lation. This behavior can be understood as following: 
For T < 60 K, the Coulomb HF term has important ef- 
fect on spin relaxation. The HF effective magnetic field 
elongates the SRT. According to Eq. (P7)) . this effect in- 
creases with the momentum scattering time. Without the 
electron-hole Coulomb scattering the momentum scatter- 
ing time is elongated, which enhances the effect and leads 
to longer SRT. For higher temperature (T > 60 K), the 
effect of the Coulomb HF term is weak)^ and the sys- 
tem returns back to the motional narrowing regime. The 
SRT thus decreases when the electron-hole Coulomb scat- 
tering is removed. The results indicate that the electron- 
hole Coulomb scattering is comparable with the electron- 
electron and electron-LO-phonon scatterings. In other 
words, besides the screening from holes, the main con- 
tribution of the hole system to electron spin relaxation 
comes from the electron-hole Coulomb scattering in in- 
trinsic semiconductors. 

In Fig. [71 we plot the SRT as function of temperature 
for iVex = 2 X 10" cm-3 with P = 2 %. In Sec. HI B2, 
we showed that there is a peak in the temperature de- 
pendence of SRT due to the electron-electron Coulomb 
scattering when the impurity density is low. In intrinsic 



B. Density dependence 

We plot the density dependence of the SRT in Fig.|5]for 
both low temperature (T = 40 K) and room temperature 
(T — 300 K) cases. It is seen that for both cases the BAP 
mechanism is far less efficient than the DP mechanism. 
Another remarkable feature is that the SRT shows a non- 
monotonic photo-excitation density dependence with a 
peak at some density n,. which resembles that in n-type 
materials. Further calculation gives Uc ~ 0.8 x 10^^ cm^"^ 
(Tf = 25 K) for T = 40 K case and = 0.9 x lO" cm"^ 
(Tf = 127 K) for T = 300 K case. Interestingly, two re- 
cent experiments give different photo-excitation density 
dependences of SRT at room temperature: in Ref. |3J 
the SRT decreases with TVex where TVex > 10"'^^ cm""^, 
while in Ref. [S^ the SRT increases with A^ex where the 
photo-excitation density is lower. These observations are 
consistent with our results. However, the peak has not 
been reported in the literature. 

We then discuss the effects of the electron-hole 
Coulomb scattering and the Coulomb HF term on the 
spin relaxation for T = 40 K as function of photo- 
excitation density. In Fig. jSj^c), we plot the SRTs ob- 
tained from the full calculation, from the calculation 
without the electron-hole Coulomb scattering, and from 
the calculation without the Coulomb HF term. It is seen 
that the Coulomb HF term plays a visible role only for 
high densities, as the HF effective magnetic field increases 
with electron density.^^ Similar to the temperature de- 
pendence [Fig. |6][c)], without the electron-hole Coulomb 
scattering, the SRT is larger for Ncx > 3 x 10"^^ cm~^ 
where the Coulomb HF term plays a prominent role, 
while it is smaller for lower photo-excitation densities 
where the Coulomb HF term is unimportant. Notably, 
the peak of SRT still exists when the Coulomb HF term 
is removed, which implies that the degree of initial spin 
polarization is irrelevant for the existence of the peak. 
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FIG. 8: (Color online) Intrinsic GaAs. SRT r due to the BAP 
mechanism and that due to the DP mechanism as function of 
photo-excitation density A'ox at: T = 40 K (a) and T = 
300 K (b) (note that the value of tbap in figure (b) has been 
rescaled by a factor of 0.01). (c): the SRT due to the DP 
mechanism from the full calculation (curve with •), from the 
calculation without the Coulomb HF term (curve with □), 
and from the calculation without the electron-hole Coulomb 
scattering (curve with A) for T = 40 K. 



V. ELECTRON SPIN RELAXATION IN p-TYPE 
III-V SEMICONDUCTORS 



In this section, we study spin relaxation in p-type III- 
V semiconductors. The main sources of spin relaxation 
have been recognized as the BAP mechanism and the DP 
mechanism;22, We first compare the relative efficiency of 
the two mechanisms for various hole densities and tem- 
peratures. After that, the hole density and the photo- 



excitation density dependences of the SRT at given tem- 
perature are also discussed. 



Comparison of the DP and BAP mechanisms in 
GaAs 



We first address the relative importance of the BAP 
and DP mechanisms for various hole densities and tem- 
peratures in GaAs. The electrons are created by photo- 
excitation with P = 50 % (i.e., we assume ideal optical 
orientation by circularly polarized light) i2i In order to 
avoid exaggerating the DP mechanism, we use the SOC 
parameter fitted from the experimental data in Ref. 0, 
i.e., 7d = 8.2 eV-A'^ (see Appendix B) throughout this 
section, which is smaller than the value from the k • p 
calculation 7d = 23.9 eV-k^M^ 



1. Low photo- excitation 

We first concentrate on low photo-excitation density 
regime, where we choose A'ox — 10^** cm^'^. The ra- 
tio of the SRT due to the BAP mechanism to that due 
to the DP mechanism is plotted in Fig. [U^a) for various 
hole densities. It is seen that the DP mechanism domi- 
nates at high temperature, whereas the BAP mechanism 
dominates at low temperature, which is consistent with 
previous investigations i^'^'^^i'^^i'*^ An interesting feature 
is that the ratio first decreases rapidly, then slowly, and 
then again rapidly with decreasing temperature. A typ- 
ical case is shown in Fig. ^ih) for = 3 x 10^* cm""^. 
It is noted that the "plateau" is around the hole Fermi 
temperature Tp = 156 K which is given by 



2kB'mo 



(71-272) ^/' + (71 + 272) 



-3/2 



2/3 ■ 



(29) 

The underlying physics is that: on one hand, the Pauli 
blocking of holes becomes important when T < Tp, 
which slows down the BAP spin relaxation effectively [see 
Eq. (pO)) ] ; on the other hand, the increase of the screening 
(mainly from holes) with decreasing temperature weak- 
ens the electron-impurity and carrier-carrier scatterings 
and thus enhances the DP spin relaxation. Consequently, 
the decrease of the ratio with decreasing temperature 
slows down and the "plateau" is formed around Tp. How- 
ever, after the hole system enters the degenerate regime, 
the screening changes little with temperature. The ratio 
thus decreases rapidly with decreasing temperature again 
[see Appendix C]. We also plot the SRT due to the BAP 
mechanism without the Pauli blocking of holes as dotted 
curve in Fig.l^I^b), which indicates that the Pauli blocking 
of holes effectively suppresses the BAP spin relaxation at 
low temperature (T < Tp). It is also seen from Fig. [H^b) 
that the total SRT increases with decreasing temperature 
as both tdp and tbap do. 
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FIG. 9: (Color online) p-GaAs. Ratio of the SRT due 
to the BAP mechanism to that due to the DP mechanism 
as function of temperature for various hole densities with 
A''cx = lO" cm~^ and rn = n,j. (a): = 3 x 10^^ cm"^ 
(curve with •), 3 x 10^^ cm~^ (curve with □), 3 x 10^'^ cm~^ 
(curve with A), and 3 x 10^* cm"^ (curve with V). The hole 
Fermi temperatures for these densities are Tp = 1.6, 7.3, 34, 
and 156 K, respectively. The electron Fermi temperature is 
Tf = 1.4 K. (b): The SRTs due to the BAP and DP mecha- 
nisms, the total SRT, together with the ratio tbap/tdp (curve 
with A) versus the temperature for Uh = S x 10^* cm~^. 
The dotted curve represents the SRT due to the BAP mech- 
anism without the Pauli blocking of holes. Note the scale of 
t"bap/t"dp is on the right hand side of the frame. 



FIG. 10: (Color online) p-GaAs. Ratio of the SRT due to the 
BAP mechanism to that due to the DP mechanism as function 
of temperature for various hole densities with A'^ex = O.luh 
and Hi — rih- (a): = 3 x 10^^ cm~^ (curve with •), 
3 X lO^'' cm"^ (curve with □), 3 x lO^'^ cm"^ (curve with A), 



and 3 x 10^* cm 



(curve with V) 

nh 



The hole Fermi temper- 
atures for these densities are Tp = 1.7, 7.7, 36, and 167 K, 
respectively. The electron Fermi temperatures are Tp = 2.8, 
13, 61, and 283 K, respectively, (b): The SRTs due to the 
BAP and DP mechanisms, the total SRT, together with the 
ratio tbap /tdp (curve with A) versus the temperature for 
rih = 3 X 10^* cm"''. The dotted (dashed) curve represents 
the SRT due to the BAP mechanism without the Pauli block- 
ing of electrons (holes). Note the scale of tbap /tdp is on the 
right hand side of the frame. 



2. High photo- excitation 



We then discuss the case with high photo-excitation 
density, where we choose A^ex = O.lrih- The ratio of 
the SRT due to the BAP mechanism to that due to the 
DP mechanism is plotted in Fig. [TOT a) for various hole 
densities. It is seen that, interestingly, the ratio is non- 
monotonic and has a minimum roughly around the Fermi 
temperature of electrons, T ^ Tp. The BAP mechanism 
is comparable with the DP mechanism only in the mod- 
erate temperature regime roughly around Tp, whereas 
for higher or lower temperature it becomes unimportant. 
To explore the underlying physics, we plot the SRTs 
due to the BAP and DP mechanisms in Fig. [TOT b) for 
n/i = 3 X 10^^ cm^''. It is seen that the Pauli blocking of 
electrons and holes largely suppresses the BAP spin re- 
laxation in the low temperature regime and hence makes 



Tbap always increase with decreasing temperature. On 
the other hand, tdp first increases with decreasing tem- 
perature, then saturates at low temperature (T < Tp) 
due to the fact that both the inhomogeneous broadening 
and the momentum scattering change little in the de- 
generate regime. Therefore, the ratio tbap/tdp first de- 
creases then increases with decreasing temperature and 
shows a minimum roughly around the electron Fermi 
temperature. This scenario holds for arbitrary excitation 
density, and the temperature where the ratio tbap /tdp 
reaches its minimum increases with excitation density. 
Finally, it is seen that the total SRT saturates at low 
temperature as tdp does. 
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FIG. 11: (Color online) p-GaAs. SRTs r due to the BAP 
and DP mechanisms together with the total SRT versus the 
photo-excitation density A'ex- The ratio of the two is plotted 
as dashed curve (note that the scale is on the right hand 



side of the frame), (a): n; = nu - 
■ 3 X 10^* cm-^ T = 50 K. 
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mechanism, it decreases slowly with the photo-excitation 
density in the low density regime, but rapidly in the high 
density regime. The decrease is mainly due to the in- 
crease of the averaged electron velocity {vk) [see Eq. (PH]) ]. 
which is determined by the temperature and is insen- 
sitive to density in the non-degenerate regime, but in- 
creases rapidly in the degenerate regime. However, the 
increase of the spin relaxation due to the BAP mechanism 
is slower than that due to the DP mechanism, because the 
inhomogeneous broadening increases as cx iV^^ while {vk) 

1 /3 

increases as oc A'ex . Consequently, the BAP mechanism 
becomes even less important in the high photo-excitation 
density regime. Similar situation also happens for other 
hole densities and temperatures. In Fig. Illf b). we plot 
the case for a larger hole density rih — iy. 10^^ cm"-^. It 
is seen that under low photo-excitation, the BAP mech- 
anism is more important than the DP mechanism. How- 
ever, the BAP mechanism becomes less important than 
the DP mechanism in the high photo-excitation density 
regime. The crossover of the low photo-excitation den- 
sity regime to the high photo-excitation density regime 
takes place around Tp ^ T. This leads to the conclu- 
sion that the BAP mechanism is not important at high 
photo-excitation density in p-type materials. It is seen 
from Fig. [Tl] that the total SRT decreases with photo- 
excitation density as both tdp and tbap do. This behav- 
ior is also consistent with what observed in experiments 
in Ref. M. 



C. Hole density dependence 



B. Photo-excitation density dependence 

We now turn to the photo-excitation density depen- 
dence of the SRT. In Fig. [Tlla), we plot the SRT due 
to the DP mechanism and that due to the BAP mech- 
anism as function of photo-excitation density for Uh = 
3 X 10" cm-3 with T = 50 K. It is seen that the SRT 
due to the DP mechanism decreases with the photo- 
excitation density monotonically. Specifically, it first de- 
creases slowly, then (Ai'ex > 10^^ cm~^) rapidly with the 
photo-excitation density. The scenario is as follows: In 
p-type semiconductors at low temperature, the dominant 
scattering mechanisms are the electron- hole and electron- 
impurity scatterings. The momentum scattering due to 
these two mechanism changes little with photo-excitation 
(electron) density (iVcx = n,,) for A'ex < n-h- In the low 
density regime {rie < 3x 10^^ cm~^, or Tp < 13 K), where 
the electron system is non-degenerate, the increase of 
density affects the inhomogeneous broadening very little. 
Thus the SRT changes slowly with the photo-excitation 
density. In the high density regime (rig > 3 x 10^^ cm^'^, 
or Tp > 61 K), the electron system is degenerate where 
the inhomogeneous broadening increases fast with den- 
sity. Consequently the SRT decreases rapidly with the 
photo-excitation density. For the SRT due to the BAP 



We also study the hole density dependence of spin 
relaxation due to the BAP and DP mechanisms. In 
Fig. [IHa), we plot the SRT due to the BAP mecha- 
nism and that due to the DP mechanism as function of 
hole density for T = 60 K and iVgx = lO^'' cm^^ It 
is seen that the SRT due to the BAP mechanism de- 
creases as l/ufi at low hole density, which is consistent 
with Eq. pop , i.e., for non-degenerate holes tbap oc l/uh- 
At high hole density, tbap decreases slower than l/rih 
due to the Pauli blocking of holes. However, for the SRT 
due to the DP mechanism, the dependence is not so ob- 
vious: the SRT first increases, then decreases and again 
increases with the hole density. As the electron distribu- 
tion, and hence the inhomogeneous broadening, does not 
change with the hole density, the variation of the SRT 
due to the DP mechanism solely comes from the mo- 
mentum scattering (mainly from the electron-impurity 
scattering). To elucidate the underlying physics, we 
plot the SRT due to the DP mechanism calculated with 
the RPA screening together with those calculated with 
the Thomas-Fermi (TF) screening^^ [which applies in 
the degenerate (high density) regime], the Debye-Huckle 
(DH) screening-^S, [which applies in the non-degenerate 
(low density) regime] in Fig. [T^ b). From the figure it 
is seen that the first increase and the decrease is con- 
nected with the DH screening, whereas the second in- 
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FIG. 12: (Color online) p-GaAs. (a): SRTs r due to the BAP 
and DP mechanisms together with the total SRT against hole 
density n^. iVex = 10^"* cm"^, rn ^ uh, and T = 60 K. The 
dotted curve denotes a fitting of the curve with • using 1 /uh 
scale. The curve with A denotes the ratio tbap/tdp (note 
that the scale is on the right hand side of the frame), (b): 
SRTs due to the DP mechanism with the DH (curve with □), 
TF (curve with A), and the RPA (curve with •) screenings. 
The ratio K^/(g^) is plotted as curve with V (note that the 
scale is on the right hand side of the frame). 



crease is connected with the TF screening. The underly- 
ing physics is as follows: In the low hole density regime, 
the screening from the holes is small and the Coulomb 
potential, which is proportional to 1/(k^ + q^), changes 
slowly with the screening constant k. Hence the electron- 
impurity scattering increases with Uh as it is propor- 
tional to ni{Vq) oc Uh (as n.i — rih). For higher hole 
density {rih > 10^^ cm"'^), the screening constant k be- 
comes larger than the transfered momentum q. [To elu- 
cidate the relative ratio of the two, we plot the ratio of 
the average of the square of the transfered momentum 
(q^) to the square of the screening constant as curve 
with V in Fig. [T2r b).] Hence the electron-impurity scat- 
tering decreases with Uh because it is proportional to 
ni{V^) oc rihl oc n~^^ as k} oc n/i for the DH screening. 
As the hole density increases, the hole system enters into 
the degenerate regime, where the TF screening applies 



and 



oc n 



1/3 



Hence, the electron-impurity scattering 



increases with the hole density as ni(y^) oc n^^' 



1/3 



Con- 



sequently, the SRT first increases, then decreases and 
again increases with the hole density as the momentum 



scattering does. It should be mentioned that this behav- 
ior is different from that in the p-type (001) quantum 
wells where tdp increases with monotonically^l as the 
screening from holes is much weaker in that case due to 
lower-dimension in phase-space and smaller hole effective 
mass [in (001) GaAs quantum wells, the in-plane effective 
mass of the heavy-hole is ~ 0.1 Imp compared to 0.54mo 
in bulk]. 

It is also noted in Fig. [T^a) that the ratio tbap / tdp 
first decreases rapidly, then slowly and again rapidly with 
the hole density n/j. The first decrease is because that 
Tbap decreases with n^, whereas tdp increases with it. 
In the crossover regime (nt ~ 10^* cm~'^), where Tp ~ T, 
the SRT due to the DP mechanism varies slowly with hole 
density. As the SRT due to the BAP mechanism also 
varies slowly with hole density due to the Pauli block- 
ing of holes, the ratio tbap/''dp changes slowly with 
hole density in this regime and a "plateau" is formed 
at Tp T. In higher hole density regime, however, tdp 
increases with n^, whereas tbap decreases with it. The 
ratio tbap/t'dp hence decreases rapidly with Uh again 
and the BAP mechanism becomes more and more impor- 
tant. It is seen from Fig. [TSJa) that the total SRT first 
increases then decreases in the low hole density regime as 
Tdp does. Consequently, the hole density dependence of 
the SRT exhibits a peak which has never been reported. 
Nevertheless, in the regime of higher hole density, the 
total SRT changes slowly with Uh as the BAP and DP 
mechanisms compete with each other. 



D. Other III-V semiconductors 

Although the above conclusions are obtained from 
GaAs, they also hold for other III-V semiconductors. 
To demonstrate that, we also investigate the problem in 
GaSb. GaSb is a narrow band gap III-V semiconductors 
of which the values of Ai^LT and AE^n can be found in 
literature In Fig.[T31 we plot the ratio of the SRT due 
to the BAP mechanism to that due to the DP mechanism 
as function of temperature for various hole densities. It 
is seen from the figures that the features are similar to 
those in GaAs, whereas the ratio is much larger than that 
in GaAs under the same condition. That is, the relative 
importance of the BAP mechanism in GaSb is smaller 
than that in GaAs. This is because the SOC in GaSb 
is much larger than that in GaAs while the longitudinal- 
transversal splitting AE'lt in GaSb is smaller than that 
in GaAs. 



VI. EFFECTS OF ELECTRIC FIELD ON SPIN 
RELAXATION IN n-TYPE III-V 
SEMICONDUCTORS 

In this section, we study the effects of electric field on 
spin relaxation in n-type III-V semiconductors. Previous 
studies have demonstrated that in quantum wells a rela- 
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FIG. 13; (Color online) p-GaSb. Ratio of the SRT due to 
the BAP mechanism to that due to the DP mechanism as 
function of temperature for ~ S x 10^^ cm~'^ (curve with 
•), 3 X 10^*^ cm"^ (curve with □), 3 x 10^^ cm"^ (curve with 
A), and 3 x 10^* cm~^ (curve with V). Ui — nn- (a): A'cx — 
10" cm-^ (b): TVcx = Q.lnu. 



lively high in-plane electric field can effectively manipu- 
late the SRT.2ii2^iS2i2a The underlying physics is that the 
high electric field induces two effects: the drift of the elec- 
tron ensemble which enhances the inhomogeneous broad- 
ening (as electrons distribute on larger k states where the 
SOC is larger), as well as the hot-electron effect which 
enhances the momentum scattering. The former tends 
to suppress while the latter tends to enhance the SRT. 
Thus the SRT has nonmonotonic electric field depen- 
dence: it first increases due to the hot-electron effect 
then decreases due to the enhancement of inhomogeneous 
broadening. In bulk semiconductors, the electric field de- 
pendence of spin lifetime has not been investigated. In 
this section, we present such a study. Using n-type GaAs 
as an example, we demonstrate that the electric field de- 
pendence of spin lifetime can be nonmonotonic (first in- 
creasing then decreasing) or monotonic (decreasing) de- 
pending on the lattice temperature and the densities of 
impurities and electrons. The underlying physics is ana- 
lyzed. The study indicates that the spin lifetime can be 
effectively controlled by electric field. 

In Fig. [141 we plot the ratio of the SRT under electric 
field to the electric-field-free one as function of electric 
field for = 2 X 10^^ cm'^ with P = 2 % at T = 150 K 
and rii = Ue- The electric field is chosen to be along 



FIG. 14: (Color online) n-GaAs. Ratio of the SRT under elec- 
tric field to the electric-field free one t{E)/t{E = 0) (curve 
with •) and the ratio of the hot-electron temperature to the 
lattice temperature T^/T (curve with A) as function of elec- 
tric field for m ^ = 2 x lO" cm"^ at T =150 K. The 
mobility is also plotted as curve with □ (note that the scale 
is on the right hand side of the frame). 



the X-axis and the initial spin polarization is along the 
2;-axis. Due to the cubic form of the SOC, the average 
of the spin-orbit field is negligible even in the presence of 
finite center-of-mass drift velocity, which is different from 
the case in quantum wells where the linear k-term gives 
a large effective magnetic field in the presence of electric 
field thanks to the strong well confinement , i^^ It is seen 
that the ratio first increases a little and than decreases 
rapidly with the electric field. At £^ = 1 kV/cm, the ratio 
drops to 0.6 which demonstrates that the electric field can 
effectively affect the SRT. To understand these effects, we 
also plot the hot-electron temperature in the figure. It is 
noted that the electron temperature increases first slowly 
then {E > 0.3 kV/cm) rapidly with the electric field, in- 
dicating clearly the hot-electron effect. For the 2DES, 
where the SOC is dominated by the linear term, the 
hot-electron effect mainly leads to the enhancement of 
scattering, whereas the enhancement of inhomogeneous 
broadening due to the hot-electron effect is marginal/^ 
Differently, the hot-electron effect also has important ef- 
fect on inhomogeneous broadening in bulk system, as the 
SOC is cubic k-dependent. As both the drift effect and 
the hot-electron effect increase the inhomogeneous broad- 
ening, the enhancement of the inhomogeneous broaden- 
ing is faster than the increase of momentum scattering. 
Consequently, the SRT is easier to decrease with the elec- 
tric field, which is different from the case of 2DES. We 
speculate that the electric field dependence of the SRT 
in Wurtzite semiconductors and strained III-V semicon- 
ductors resembles that in the GaAs quantum wells, when 
the SOC is dominated by the linear k termJ^ii^iSaiSa The 
mobility of the electron system is also plotted in the fig- 
ure. The variation of the mobility indicates the nonlinear 
effects of the electric field in the kinetics of electron sys- 
tem. 

For lower and higher temperature cases, we plot the re- 
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FIG. 15: (Color online) 7i-GaAs. Ratio of the SRT under elec- 
tric field to the electric-field-free one t(E)/t{E = 0) (solid 
curves) and the ratio of the hot-electron temperature to the 
lattice temperature T^/T (dashed curves) as function of elec- 
tric field for (a): T = 30 K with rie = 10^'' cm"^ (curve with 
• : Ui = rie; curve with □: Ui — 0.05ne) and (b): T = 300 K 
with He = rii (curve with •: rie = 10^^ cm~'^; curve with □: 



suits in Fig. [151 The electron density is Ue = 10^® cm^"^ 
in Fig.fTSTa). It is seen that the ratio for the case with low 
impurity density rii — O.OSrte decreases faster than that 
for the case with Ui — rie- This is because both the hot- 
electron effect and the drift effect are more pronounced 
in cleaner system,^** which thus leads to a faster decrease 
of the SRT due to the enhancement of inhomogeneous 
broadening. In the background of hot-electron spin in- 
jection under high bias, our results indicate that the ma- 
nipulation of the SRT by the electric field is more pro- 
nounced for samples with high mobility. For high temper- 
ature case {T — 300 K), we plot the SRT and hot-electron 
temperature as function of electric field for two different 
electron densities Ue = 10"'^^ cm~^ and rig = 2x 10^"^ cm^'^ 
with Tii = Tie in Fig. [T5f b). It is seen that both the SRT 
and the electron temperature differ marginally for the 
two cases even though their impurity and electron densi- 
ties differ by 20 times. This is because that at 300 K the 
electron-LO-phonon scattering is more important than 
the electron-impurity scattering. Thus both the drift of 
the electron system and the hot-electron effect is mainly 
determined by the electron-LO-phonon scattering, and 
the variation of the SRT with electric field is insensitive 



to impurity density. For electron density, as the electron 
system is in the non-degenerate regime for both cases, 
the electron density dependence is hence also weak. 



VII. CONCLUSION 

In conclusion, we have applied an efficient scheme, the 
fully microscopic KSBEs, to study the spin dynamics in 
bulk III-V semiconductors, with all scatterings explicitly 
included. This approach takes full account of the spin 
relaxation due to the DP, EY, and BAP mechanisms in 
a fully microscopic fashion, and enables us to find impor- 
tant results missing in the previous simplified approaches 
in the literature. From the KSBE approach, we study the 
electron spin relaxation in n-type, intrinsic, and p-type 
III-V semiconductors. We also investigate the effects of 
electric field on spin relaxation in rt-typc III-V semicon- 
ductors. 

For n-type III-V semiconductors, differing from the 
previous conclusions, we find that the spin relaxation due 
to the EY mechanism is less important than that due to 
the DP mechanism even in narrow band-gap semicon- 
ductors, such as InAs and InSb. We then focus on the 
spin relaxation due to the DP mechanism. We find that 
the density dependence of the SRT is nonmonotonic and 
we predict a peak which appears in the metallic regime. 
This behavior is due to the different density dependences 
of the inhomogeneous broadening and the momentum 
scattering in the degenerate and non-degenerate regimes. 
In the non-degenerate regime, as the electron distribu- 
tion is close to the Boltzmann distribution, the inhomo- 
geneous broadening changes little with the density but 
the electron-electron and electron-impurity scatterings 
increase with the electron density. As a result, the SRT 
increases with the density. In the degenerate regime, the 
inhomogeneous broadening increases with electron den- 
sity, whereas the momentum scatterings decrease with it. 
Consequently, the SRT decreases with the electron den- 
sity in the degenerate regime. The peak of the SRT is 
hence formed in the crossover regime, where the corre- 
sponding Fermi temperature is close to the lattice tem- 
perature, Tp ~ T. Our results show that the electron- 
electron scattering plays an important role in the spin 
relaxation. We also study the density dependence for the 
case with strain-induced SOC, where the density depen- 
dence of the inhomogeneous broadening is different due 
to the linear k-dependence of the strain-induced SOC. 
However, the SRT still has a peak but at a larger den- 
sity compared to the strain-free case. We further study 
the temperature dependence of the SRT. We find that 
the SRT decreases monotonically with the temperature 
which is consistent with experimental findings. After we 
artificially lower the impurity density, we find a peak in 
the SRT which is due to the different temperature depen- 
dence of the electron-electron scattering in the degener- 
ate and non-degenerate regimes. This is consistent with 
the results in 2DES where the peak in the SRT due to the 
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electron-electron scattering appears only when the impu- 
rity density is low (e.g., = O.lnf.)^ We also study the 
initial spin polarization dependence of the SRT where 
the effect of the Coulomb HF term is discussed. We find 
that the dependence is quite weak in bulk system com- 
pared to that in the 2DES, which is again due to the 
large impurity density rii > Ue in bulk system. 

For intrinsic III-V semiconductors, we first compare 
the BAP mechanism and the DP mechanism. We find 
that the BAP mechanism is far less efficient than the 
DP mechanism. We further compare our results from 
the fully microscopic KSBE approach with those from 
the approach widely used in the literature. We find 
that the previous approach deviates in the low tem- 
perature regime due to the pretermission of the Pauli 
blocking. Also, the previous approach ignores the long- 
range electron-hole exchange scattering which is shown 
to be dominant in GaAs. We find that the electron-hole 
Coulomb scattering plays an important role in spin relax- 
ation. The Coulomb HF term is found to have important 
effects on spin relaxation at low temperature and high 
photo-excitation density, as the impurity density is very 
low in intrinsic semiconductors (we choose rii = 0). Due 
to the same reason, the peak in the temperature depen- 
dence of the SRT due to the electron-electron scattering 
also appears at small spin polarization. We further dis- 
cuss the photo-excitation density dependence of the SRT. 
We find that the SRT first increases then decreases with 
the density which resembles the case in n-type samples 
where the underlying physics is also similar. 

For p-type III-V semiconductors, we first examine the 
relative importance of the BAP mechanism. We find that 
the BAP mechanism dominates the spin relaxation in the 
low temperature regime only when the photo-excitation 
density is low. However, when the photo-excitation den- 
sity is high, the BAP mechanism can be comparable 
with the DP mechanism only in the moderate temper- 
ature regime roughly around the Fermi temperature of 
electrons, whereas for higher or lower temperature it is 
unimportant. The photo-excitation density dependences 
of SRTs due to the BAP and DP mechanisms are also dis- 
cussed. We find that the relative importance of the BAP 
mechanism decreases with photo-excitation density and 
eventually becomes negligible at sufficiently high photo- 
excitation density. For hole density dependence at small 
photo-excitation density, we find that the spin relaxation 
due to the BAP mechanism increases with hole density 
linearly in low hole density regime but the increase be- 
comes slower in high hole density regime where the Pauli 
blocking of holes becomes important. Interestingly, the 
SRT due to the DP mechanism first increases, then de- 
creases and again increases with hole density. The under- 
lying physics is that the momentum scattering (mainly 
from the electron-impurity scattering) first increases with 
hole (impurity) {rii — n^) density, then decreases with 
hole density due to the increase of the screening. How- 
ever, at high hole density when the hole system is degen- 
erate, the screening increases slower with the hole density 



and the momentum scattering again increases with the 
hole (impurity) density. On the other hand, the inhomo- 
geneous broadening does not change with hole density as 
the electron density is solely determined by the photo- 
excitation density. Consequently, the SRT due to the 
DP mechanism first increases, then decreases and again 
increases with the hole density. This behavior makes the 
ratio tbap /^dp first decreases rapidly, then slowly and 
again rapidly with the hole density. The BAP mecha- 
nism is more important than the DP one for high hole 
density. The relative importance of the BAP mechanism 
in GaSb is found to be much less than that in GaAs due to 
both the weaker electron-hole exchange interaction and 
the larger SOC in GaSb. 

Finally, we study the effect of electric field on the spin 
relaxation in n-type GaAs. We find that the SRT can 
be largely affected by the electric field. The underlying 
physics is that the electric field induces two effects: the 
center-of-mass drift which enhances the inhomogeneous 
broadening and the hot-electron effect which increases 
both the momentum scattering and the inhomogeneous 
broadening. The electric field dependence of SRT thus 
can be nonmonotonic: it first increases due to the in- 
crease of scattering then decreases due to the enhance- 
ment of inhomogeneous broadening. However, we find 
that differing from the 2DES, the SRT is easier to de- 
crease with the electric field. This is because that the in- 
homogeneous broadening increases faster when the SOC 
is cubic compared to the 2DES where the SOC is domi- 
nated by linear k term. We expect that the electric field 
dependence of the SRT resembles the 2DES for Wurtzite 
semiconductors or strained semiconductors, where the 
SOC can be dominated by the linear k term. We also 
find that the effect of the electric field becomes more 
significant for low impurity density samples at low tem- 
perature as both the drift effect and the hot-electron ef- 
fect are more pronounced. However, at room tempera- 
ture, the effect of the electric field is insensitive to the 
impurity density as the electron-LO-phonon scattering 
is more important than the electron-impurity scattering. 
The electron density dependence is also weak as long as 
the system is in the non-degenerate regime. 

Note added: After we submitted this manuscript, the 
peak we predicted in the doping density dependence of 
the SRT of n-GaAs in the metallic regime has been re- 
alized experimentally in a subsequent paperiiSS, Also an 
independent theoretical calculation from the KSBE ap- 
proach well reproduced the peak in the same paper 
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APPENDIX A: NUMERICAL SCHEME 

Our numerical scheme is based on the discretization of 
the k-space similar to that in Ref. 113. Here we extend 
it to the three dimensional case. Our technique greatly 
reduces the calculation complexity and makes the quan- 
titatively accurate calculation possible. 

The k-space is divided into N x M x L control re- 
gions where the k-grid points are chosen to be \in.m,i = 
\/2mcEn(sm 9m cos 0; , sin 9^ sin , cos 9m)- To facilitate 
the evaluation of the (5- functions in the scattering terms, 
we set En — (n + 1/2)AE, where the energy span in 
each control region is AE = wlo/'t-lo with tt-lo being 
an integer number and wlo denoting the LO-phonon fre- 
quency. The electron-impurity and electron-phonon scat- 
terings are then solved easily since the (5-functions can be 
integrated out directly. 

For the electron-electron Coulomb, electron-hole 
Coulomb and electron-hole exchange scatterings, the sit- 
uation is much more complex. The electron-electron 
Coulomb scattering term [Eq. ([9])] can be rewritten as, 

^^P^\ee - E(2;^^k-k'[Ak,k'P^,Ak^k/5^i^(k,k') 

- AkM'Pk'^k'MP^H{k\ k)] + H.C., (Al) 



with 9o = arccos[(g^ -I- uj)/2k"q]. Note that the inte- 
gration over Ek" is restrained by the condition £k" ^ 
[(g^ -f a;)/2q]^/2TOc according to |cos6'op < 1. Now the 
electron-electron Coulomb scattering is easily integrated 
out. Note that the calculation of H{q,uj) can be done 
before the calculation of the electron-electron Coulomb 
scattering terms, which thus reduces the whole calcula- 
tion complexity from 0{N^APL^) to 0{N^M^L'^). This 
method, first developed by Cheng in two-dimensional 
system,'i2i greatly reduces the calculation complexity. 

For the electron-hole Coulomb scattering, the idea is 
similar, but the technique is more complex. Denoting 
Xm = "'nc/nim, the (5- function in Eq. (|10p can be written 
as 

S{ek"m - ek"-qm' " w/2mc) = 2mc 
x<^ (^{xm — Xm') k"^ — Xm'Q^ + 2xm'k" qcos9 — iSj .(A6) 

The (5-function is then integrated out similarly. 



(new) 



(A7) 



The integration over is restrained by the condition 
I cos6'oP < 1, where 



where 



1 



i?(k, k') = (2^)3 ^ Siek" - £k"-k+k' + ek' - £k) 



xTr (^Ak",k"-k+k'/5k"_k+k'Ak"-k+k'.k"/5k" j ■ (A2) 
Substituting q = k — k' and uj = 2mc(ek — Ek'); '^^^ 
iJ(q,cj) y dk"(5(ek" -£k"-q-t^/2mc) 



XTr(^Ak",k"-qP^,,_qAk"-q,k"/5^.J . (A3) 

Now the 6 function can be simplified as 
S{s-k" - ek"-q - w/2mc) = —^S{cos9- cos 6*0), (A4) 

where 9 is the angle between k" and q, and cos 6*0 = 
(g^ +uj)/{2k"q). To evaluate the 5 function, it is helpful 
to rotate to the new coordinate system with q being along 
the z-axis. In this coordinate system, 9 — 9 and the S 
function can be evaluated readily. The result is 

H{q,Lu) = / d£k"rf0" 

q J 



XTr( Ak",k"-q/5k-_qAk"-q,k"/5k'. 



'new) 



,(A5) 



COSfo = 



2k" q \x 



(A8) 



The restriction can be simplified as 



b^k"^ -{1 + 2ab)k" < and k" > 0, (A9) 

where a = +q^)/{2q) and b = (a;„ - Xm') / {2xm'q)- 
The above inequality can be worked out readily for given 
{q,uj,m,m'), and the restriction condition for = 

Xmk""^ / {2mc) is then obtained. 

The electron-hole exchange scattering is solved simi- 
larly by substituting the matrix element I'^k'/Z^^'P ^'^^^ 

Finally, the drift term is solved with similar 
method of that in Ref. [13. The differential equations are 
solved by the fourth-order Runge-Kutta method. 

The computation is carried out in a parallel manner by 
using OpenMP. For a typical calculation with the parti- 
tion of 40 X 8 X 16 grid points in the k-space, it takes 
about ten hours to evolute 50 ps on a quad-core AMD 
Phcnom 9750. 



APPENDIX B: COMPARISON WITH 
EXPERIMENT 

We compare the calculation from the fully microscopic 
KSBE approach with the experimental results in Refs. [3 
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FIG. 16: (Color online) (a) n-GaAs. SRTs r from the exper- 
iment in Ref. (•) and from the calculation via the KSBE 
approach with only the DP mechanism (solid curve) as well 
as that with only the EY mechanism (dashed curve), — 
W^^ cm"^ m = rie and iVox = 10^" cm"^ 70 = 8.2 eV-A^ 

(b) p-GaAs. SRTs r from the experiment in Ref. (•) and 
from the calculation via the KSBE approach (solid curve). 
Tih = 6 X 10^^ cm-^*, Hi = Uh and T = 100 K. 70 = 8.2 eV-A^. 

(c) p-GaAs. SRTs r from the experiment in Ref. (•) and 
from the calculation via the KSBE approach (solid curve). 
Uh = 1.6 X 10^*^ cm-^ m = Uh and TVex = 10^* cm'^ 
7D = 10 eV-A^ 



and [95I . These experiments were carried out in n- and 
p-type GaAs respectively. In Fig.fWa). we plot the SRT 
as function of temperature calculated from the KSBEs, 
together with the experimental data in Ref. for n-GaAs 
with Tie = 10^^ cm~^, Ui = Ue and TVex = 10^'' cm~^. It 
is seen that the calculation agrees well with the experi- 



mental results in n-GaAs for T > 20 K. The deviation in 
the lower temperature regime is due to the rising of the 
localization of electrons ii^ii^iSi The SRT due to the EY 
mechanism is also plotted in the figure, which is much 
larger than the experimental data, indicating the irrele- 
vance of the EY mechanism. The calculation gives a fit 
of the SOC parameter as 70 = 8.2 eV-A^ which is dif- 
ferent from the value 70 = 23.9 eV-A'^ calculated from 
the tight-binding or k • p parametric theoriesi^i^ How- 
ever, the value is still in the reasonable range of 70 [for 
lists of 7d calculated and measured via various meth- 
ods, see Ref. Our fitting, though with only one 
fitting parameter 7d, agrees well with the experimen- 
tal data in almost the whole temperature range, which 
is much better than the fittings with the same experi- 
mental data in Refs. [tiHtI In Fig. \Wlh), we plot the 
photo-excitation density dependence of the SRTs from 
our calculation and from the experiment in Ref. 95 for p- 
GaAs with nn = 6 x 10^^ cm'^, n, = nu and T = 100 K. 
It is seen that with the same 7d, the calculation again 
agrees well with the experimental data. The SRT due 
to the BAP mechanism is about 20 times larger, which 
is consistent with our conclusion that the BAP mecha- 
nism is unimportant at high photo-excitation density. In 
Fig. fTHT c). we plot the temperature dependence of the 
SRTs from our calculation and from the experiment in 
Ref. Ho for p-GaAs with nu = 1.6 x 10^^ cm^^, n.j = Uh 
and iVox = 10^^ cm^'^. The best fitting gives a slightly 
larger 7d = 10 eV-A^. The SRT due to the BAP mecha- 
nism is ~100 times larger than that due to the DP mech- 
anism. This is consistent with our conclusion that the 
BAP mechanism is unimportant at high temperature for 
low doping density. 

Throughout the paper, the SRT in GaAs due to the 
DP mechanism is calculated with 70 = 23.9 eV-A'^ un- 
less in Sec. V where we use 7d = 8.2 eV-A^^ in order to 
avoid possible exaggeration of the DP mechanism by us- 
ing a "larger" SOC parameter. However, the SRT due 
to the DP mechanism is proportional to 7^ in motional 
narrowing regime, thus the results presented in the pa- 
per can be easily converted to the results for another 7d. 
The ratio for the SRTs in the two cases is 8.5. 



APPENDIX C: ROLE OF SCREENING ON DP 
SPIN RELAXATION IN p-TYPE GAAS 

We study the effect of screening on DP spin relaxation 
in p-type GaAs. We focus on the temperature depen- 
dence of the SRT which corresponds to the discussions on 
the results in Figs.[9landfT0l As in Sec. V, we discuss both 
the low and high excitation density cases. We first study 
the low excitation case. To elucidate the role of screening 
on the DP spin relaxation, we plot the SRT due to the 
DP mechanism with the TF screening^^ [which applies 
in the degenerate (low temperature) regime], the DH 
screening^^ [which applies in the non-degenerate (high 
temperature) regime], and the RPA screening^^ (which 
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FIG. 17: (Color online) p-GaAs. The SRTs due to the DP 
mechanism with the TF (curve with □), DH (curve with A), 
and the RPA (curve with •) screenings for excitation density 
(a) A^ex = 10" cm"^ and (b) Aox = 3 x lO" cm'^ Hole 
density nn ~ S x 10^* cm~^ and rii — nu- 



applies in the whole temperature regime) in Fig. [TTT a). 



The crossover from the DH screening to the TF one is 
clearly seen in the figure. Also, compared with the case of 
temperature-independent screening (i.e., the TF screen- 
ing), one can see that the increase of the SRT with de- 
creasing temperature is indeed slowed down around Tp 
(where the screening effect is prominent) due to the in- 
crease of screening. 

For high excitation case, we find a peak in the tem- 
perature dependence of the SRT due to the DP mecha- 
nism. This peak is roughly around Tp. To explore the 
underlying physics, we also plot the SRT due to the DP 
mechanism with the TF and DH screenings. It is seen 
that with the TF screening the peak disappears, while 
with the DH screening the peak remains. This eluci- 
dates that the appearance of the peak is due to the in- 
crease of screening at low temperature. What is different 
for the high excitation case is that the electron Fermi 
temperature is higher than the hole Fermi temperature. 
Thus at low temperature when the electron system is 
in degenerate regime (hence inhomogeneous broadening 
changes slowly) and the hole system has not yet entered 
into the degenerate regime (hence the increase of screen- 
ing with decreasing temperature is prominent), the SRT 
due to the DP mechanism tdp decreases with decreas- 
ing temperature. On the other hand, at high tempera- 
ture, Tdp increases with decreasing temperature. Con- 
sequently, the peak is formed. It should be noted that 
this peak is different from the peak found in intrinsic 
semiconductors in Sec. IV where the peak is due to the 
nonmonotonic temperature dependence of the electron- 
electron Coulomb scattering, whereas the main scattering 
mechanism here is the electron-impurity scattering. 
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